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Abstract 

| We calculate the leading twist valence quark distribution in the pion in the framework of QCD 

! sum rules with nonlocal condensates. Particular attention has been paid to the correct account 

^ ■ for the bilocal power corrections. 
O ' 

pu,: 

> 

• r- 1 

X 



1. One of the most important achievements of Quantum Chromodynamics is the determi- 
nation of Q 2 -evolution law for the structure functions Fi(x,Q 2 ) of deep inelastic scattering. It 
allows one to calculate the magnitudes of observable Fi(x, Q 2 ) at some scale Q 2 starting from its 
value at another one Qq. The theoretical basis for this application is provided by factorization 
theorems which give a possibility to express the physical cross section as a "hard" parton sub- 
process convoluted with a "soft" parton distribution function. While the former can be treated 
perturbatively due to the celebrated asymptotical freedom, the latter is governed by strong in- 
teraction dynamics at large distances and, therefore, it cannot be evaluated within perturbative 
QCD. On the other hand, the deep inelastic cross section can also be computed by using the 
operator product expansion (OPE). This gives structure functions in terms of certain coefficients 
multiplied by the target matrix elements of local quark and gluon operators of definite twist. 
Combining the two approaches allows one to express the parton distributions in terms of quark 
and gluon correlation functions on the light cone. Following Collins and Soper we can write 
for the twist- 2 valence quark distribution in a hadron 



(h(p)\o(^n,~n \ \h(p)) = 4 jf cos(Xx)u h (x 

whereQ 
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u(~n)+(A->-A) (2) 



and $ is a path ordered exponential in the fundamental representation of the colour group along 
the straight line which insures the gauge invariance of the parton distribution 

$[x, y] = P exp (ig(x - y)„ £dat a B;(y + a(x - y))) . (3) 

It should be noted that the light-cone position representation is useful to make contact with 
the OPE approach while the light-cone fraction representation is appropriate for establishing the 
parton language. 

The determination of parton distributions is, up to now, reserved to experimental studies 

but as a final goal they are expected to be evaluated from the first principles of the theory. In 

the lack of complete understanding of the yet unclear confinement mechanism they provide a 

challenging task for nonperturbative methods presently available. Among the approaches which 

account for nonperturbative effects the most close to QCD perturbation theory are the QCD sum 

rules 0. In the last decade they were applied with moderate success to determine nucleon and 

photon structure functions in the region of intermediate values of the Bjorken variable [§, || and, 

recently, in the small A-region for the light-cone position representation [f|. 

1 Throughout the paper + subscript means the convolution of the corresponding Lorentz index with light-cone 
vector n^, such that n 2 — 0, (np) — 1, (nq) — and q is t-channel momentum introduced below. 
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While the nucleoli structure functions are now well defined by the analyses of the precise 
experimental data and are attacked theoretically, much less is known about the parton distribution 
of other hadrons, in particular of 7r-meson. Being of interest in their own right they provide good 
testing ground for predictions of the QCD sum rule method which will be used in the following 
for the determination of the leading twist pionic valence quark distribution. 

2. In order to evaluate the quark distribution in the pion by means of the QCD sum rules 
method, we consider an appropriate three-point correlation function of two axial currents that 
have non-zero projection onto the pion state being proportional to the pion decay constant 
(Oli^kG )) = ifnPfi and the nonlocal string operator O on the light cone defined by eq. (fj) 

W^( Pl ,p 2 , q)=i 2 j dWye^ x+ ^(0|T jj'J(x), o(y + ±n, y - ^ , fj (0) j |0>. (4) 

The usual strategy is to use the duality between the hadronic and partonic representation for the 
correlator under investigation. 

On the one hand, we should consider the dispersion relation for the latter and extract the 
contribution due to the low lying hadron, namely, due to 7r-meson, approximating the higher 
state contribution by perturbative spectral density 

4/ 2 f 1 

W ++ (pl,pj,Q 2 = 0) = — — — — \ dx cos(Xx)u 7T (x) 

(pi - ml)(pi - ml) Jo 

+~2 / / ds x ds 2 —^ — 2T7 2t( 1 -0(so-Si)0(s o -s 2 )), ( 5 ) 

7T 2 Jo Jo (Si - P l)(s 2 - P i) 

with parameter s characterizing the beginning of the continuum. Note that projecting the 
Lorentz indices of the pion interpolating fields on the direction picked by vector n^, we extract 
the leading tensor structure in the infinite momentum frame. We omit the subtraction polynomials 
in p 2 and p 2 because they disappear after the Borel procedure has been applied. The latter leads 
to exponential suppression of the excited state contribution in the phenomenological side of the 
sum rule and gives factorial improvement of the OPE series at the theoretical one. We perform 
the double Borel transformation and put the parameters equal M\ = Mf = 2M 2 in order not 
to introduce the asymmetry between the initial and final pion states and to make contact with 
two-point sum rules for the pion decay constant. 

On the other hand, we consider the OPE for the same quantity. Of course, the QCD sum rules 
with local condensates are inappropriate here because the usual local power corrections produce 
5-type contribution to the distribution function. It is not surprising since some propagators are 
substituted by constant factors that do not allow the momentum to flow and the whole hadron 
momentum be carried by a single quark. The probability density of this configuration in the 
phase space is 5(1 — x). Higher condensates produce even more singular terms. However, this 
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singular contribution can be smeared over the whole region of the momentum fraction from zero 
to unity by avoiding the Taylor expansion of the generic nonlocal objects which are the starting 
point of all QCD sum rule calculations and introducing the concept of nonlocal condensate |6], |7| 
which assumes the finite correlation length for the vacuum fluctuations. 

At the two-loop level, to which we restrict our analysis, we need the bilocal quark and gluon 
condensates, trilocal quark-gluon condensates and four-quark condensates. The latter will be 
factorized into the product of bilocal scalar quark condensates via the vacuum dominance hy- 
pothesis. For explicit calculations, it is convenient to parametrize the bilocal condensates in the 
form of the well-known a-representation for propagators |7j 

(0$(0)$[0,x]V(x)|0) = (#) f°° daf s (a)e ax2 /\ 

Jo 

(0$(0)$[0,x] 7m V(x)|0) = -ix M Jw^> 2 r daf v {a)e a ^/\ (6) 

oi JO 

One comment concerning eqs. (P) is that in deriving a QCD sum rule one can always perform 
a Wick rotation x — > ix and treat all the coordinates as Euclidean, x 2 < 0. We use the following 
ansatz for the distribution of vacuum quarks in the virtuality a || 

fs{a) = 2KK§KJf) exp (-7T ' ai ) ' (7) 

which gives the exponential fall-off for the coordinate dependence of the condensates found on 
a lattice 0. Here A 2 = 0.2GeV 2 , and 7 is fixed from the lowest nontrivial moment of the 
distribution function f$ that is related to the value of the average virtuality of the vacuum 
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3. Conventional calculations of the perturbative diagram (see fig. 1(a)) with the light quark 
masses neglected result in 

y pert (M 2 ,Q 2 ,\) = ^-M 2 f 1 dxcos(Xx)xxexp [ j . (8) 
Zn z Jo \ x4:M A J 

Note that we have kept the t-channel momentum transferred to be nonzero. If we expand the 
cosine in the Taylor series and integrate over x, we find out that each moment possesses logarithmic 
non-analyticities of the type (Q 2 ) n \nQ 2 . These terms come from the small-x region, where the 
spectator quark carries almost the whole momentum of the pion, so that the struck quark becomes 
wee and can propagate over large distances in the t-channel. Therefore, we have to perform 
additional factorization for separation of small and large distances in the corresponding invariant 
amplitude; this will lead to the appearance of additional terms in the OPE for the three-point 
correlation function which correct the small-x dependence of the parton density. 
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The simplest nonperturbative correction comes from the vector condensate (fig. 1(b)) 

8 

\EV(M 2 ,A) = — Tca s (uu) 2 / dx cos(Xx)x f v (xM 2 ). (9) 
81 Jo 

The dominant contribution is due to the four-quark condensate. For calculation of two loop 
diagrams appearing in the consideration (fig. 1(c)) it is very convenient to use the following 
method which is an extension of the calculation technique developed in ref. |TTJ for two-point 
correlators. The main ingredient is a construction of a more general object, namely, the current 
in the vertex opposite to the gluon propagator should be replaced by the nonlocal one with light- 
like separation. The advantage of this substitution results in appearing of extra 5-function and 
introducing through this replacement a set of variables which give the simplest integration. At 
the end, we put the nonlocallity parameter equal zero. Performing straightforward calculations 
we obtain 

32 ,.,/■! 



62 2 f 

\Ps(M , A) = — iia s {uu) / dxcos(Xx) 
9 J o 



dyl dtldCf 8 ( TM 2 )fs( |m 2 W LS) *» (10) 



A J \C J \x - v J \ny£C - xy£C\' 

The gluon as well as trilocal quark-gluon condensate contributions are numerically much less 
important than the power correction we accounted for; therefore, we neglect them in what follows. 

It is well known that there exists a parton sum rule that implies that the pion contains one 
w-quark. Summing the calculated contributions and taking the formal limit Q 2 — > in the 
perturbative term we can convince ourselves, comparing the result with the sum rule for the pion 
decay constant^, that the normalization condition is broken. The reason for this has already been 
mentioned earlier and we elaborate this point below. 

4. Now we derive a Ward identity (WI) || and show that the parton sum rule should be 
exact in QCD. Of course, from the fact that O is a point-splitted vector current it follows that in 
the limit A — > the correlator (||) is related to the derivative of the two-point correlation function 
of two axial currents. However, a more general WI (for arbitrary A) will be useful in the following 
for discussion of the condensate contribution omitted. 

Noticing that we are interested in the limit Q 2 = 0, we choose vectors and to be 
proportional. Then, integrating by parts in eq. (^j) and using the equation for the complete 
variation of the phase factor with respect to the smooth variation of the path r — > V : x M (r) — > 
xV(<r)=x M (T) + 5 X/i (T) 

5$r[x,y] = <Mx',y'] - $ r [x,y] = igt a B;(x)6x^ r [x, y] - ig$ T [x, y]*°S»(y)<Jx M (0) 

r 1 dx (t) 

+ig J^rlxMrWG^Mr^Mr)^ 1 ^^)^} (11) 

2 Of course, the comparison should be made with the sum rule by accounting for nonlocal condensates given in 
ref. §. 
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with x M (l) = x M , x M (0) = y M , we obtain 

^(pi,p 2 ,g) 



ci 4 xe ipix 



X i 



(gx) 



5t, 



+^ /V/ ^V pix+i9y ^y (0|T { #(x), S (y + ^ 



"■v - ^n,r bf(0) |0) + (A^-A). 



Here 



1 Y 2 n ' Y ~ 2 n ' T 



(12) 



A 



A 



y + 2 n '^ + r 2 n 



Tp^+fy + r^nU 



A A 



(13) 



and 



£(x,x- An) = rf(x)7^7 5 $[x,x - An]w(x - An) 



(14) 



is a point-splitted operator which when sandwiched between the pion state and that of the vacuum, 
and convoluted with the light-like vector n M defines the leading twist-2 pion wave function. From 
this WI it follows that the normalization of w-quark distribution is exact in QCD, provided it is 
not spoiled by continuum subtraction 



/ dxu 7r (x) 
Jo 



(15) 



5. As we have seen, in the limit Q 2 — > the perturbative term though finite contains the 
logarithmic non-analyticities at this point. This is a typical example of the mass singularities in 
the QCD sum rules framework |13|, |14|, 15| . In order to get rid of this perturbative behaviour and 
replace it by a physical one, it is necessary to modify the original OPE. For the form factor type 
problem a two-fold structure of the modified OPE has been realized in refs. [16|, [T^J being of the 
following schematic form: 



W(pl,pl q 2 ) = £ C^iplpl q 2 )(O d ) + £ / rf 4 xe^ x C«(x)>V l (g,x, A). 



(16) 



An additional second term determines the contribution due to the long-distance propagation of 
quarks in the t-channel. Here W are the two-point correlators 



W i (g,x,A) = J Ae w (0|Tjo i (x,0),o(V+^ 



(17) 
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of the operator in question and some nonlocal string operator of a definite twist [15 1 which arises 
from the OPE of T-product of pion interpolating fields 

T{jU*)J & J (0)} = £C W (*)^M). (18) 

i 

The coefficients C^ d \p\,p\, q 2 ) in eq. ([16]) are free from non-analyticities or singularities in Q 2 be- 
cause they are defined as the difference between the original diagram and its factorized expression 
which is the perturbative analogue of the corresponding bilocal correlator. The bilocals cannot 
be directly calculated in perturbation theory but we can write down the dispersion relation for 
them 

w( g ,M)=-r ^' H f ,A) . ( i9 ) 

7i Jo s — q l 

accepting the conventional spectral density model: "low- lying hadron plus continuum". The 
parameters of the model could be found from auxiliary sum rules. There is no need in additional 
subtractions in eq. (|l]|) because one always deals with the difference between the "exact" bilocal 
and its perturbative part; so due to the coincidence of their UV behaviours the subtraction terms 
cancel in this difference. 

6. The simplest bilocal power correction is given by the following convolution: 

W£l(phpl Q 2 , A) = |d 4 xe^ x C| 1) (x)W]: + (g,x, An), (20) 

where the coefficient function is expressed through the quark propagator CW(x) = 25 + (x) and 

WX + (q,x,Xn)=z |Ae i9y (0|T|w(0) 7+ $[0,x] M (x), O (y + ±n,y - ±n)\\Q). (21) 



We extract the contact term |[T7| , |18|1 due to the vector condensate from this correlator and saturate 
the remaining part[] by the contributions of the mesons of increasing spin; these are p°, g states 
and so on. The fact that we are interested in the C-odd distribution (valence quark) results in 
contribution of spin-odd states in the t-channel. It is very convenient to parametrize the appearing 
matrix elements via the wave functions describing the light-cone momentum fraction distribution 
of quarks inside mesons. To the leading twist accuracy we can write 

(0|^(0)$[0,x] 7 ^(x)|M J (g,r / )) = e^ 1 ^ J _x, 1 x, 2 ...x, J _ 1 (m M ) J ff ty}\xq) 

~ ^ e ^2...w x Mi x ^-- x M./ (^m) J /j 2) 0j 2) N), (22) 



3 Although for the present problem the calculation of this part is only of academic interest, as it vanishes in 
the forward limit being proportional to Q 2 , we nevertheless evaluate it in order to demonstrate the difficulties one 
faces when the contact-type contribution is absent and the estimation of the correlator is carried out saturating it 
by contribution of physical states. 
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where J is a spin of the meson, r] its polarization and e^] lM2 ^ J _ 1 is a polarization tensor. Inspired 
by our knowledge that in some cases the asymptotical wave functions turn out to be rather close 
to "exact" ones, we take in our estimation the former in the following forms that are governed by 
conformal arguments (/3 = 1 — 0) 

In our model for the bilocal correlator we can achieve this result if we assume the duality 
of the meson resonances to the bare quark loop. In general, this quite severe assumption turns 
out to be reasonable for the case at hand, at least for the mesons of the lowest spins J. It is 
known experimentally that the physical cross section averaged over the p-meson peak coincides 
with the quark one. Local duality for the low lying states is a nontrivial dynamical property and 
is not realized in all channels | l9fl . For the problem at hand it can be explained by the specific 



interaction of the classical vector mesons with the quark and gluon condensates [0. The power 
correction for them even at M 2 « m 2 does not exceed 10 — 20% of the main perturbative term. 
So, p is predicted to be dual to the quark loop with the duality interval about s p ~ 2m 2 p . However, 
the local duality for the higher spin mesons can be broken f20fl . 

The net result for the difference between the "exact" bilocal and its perturbative part reads 



W^ + (g,x,An)-W^ ert) (g,x,An) 
= zx + jf* dre ir ^ | (0|«^nj 7+ $ 

+2Q 2 (x + ) 2 / dr df3pe iT ^ qx) V ^ 2) (/3)(-Ax 
Jo Jo T j~ri> 



A A 

— n, x n 

2 2 



u(x- -n)\0) + (A -f -A) 



2 



■7=1,3,. 



87r 2 2 J - 1 r(J)r(2J + 4)y s + Q 2 v ' m 2 M + Q 2 

where <Tj is the continuum threshold and tum is the mass of the lowest meson state in the channel 
of given spin J. Substituting this expression into eq. ( p0|) and performing simple calculations we 
obtain (diagrammatical representation is given in fig. 1(d)) 

^[(M 2 ,A) = — 7ra s (uu} 2 J dx cos(Xx)x f v (xM 2 ) +Q 2 e^ £ (iA) J_1 



i7i 2 T(J)T(J + 2)Jo V ' \M 2 ) m 2 M + Q 2 J 

where Zj = <jj/4M 2 . The former term is a contact-type contribution due to the vector condensate. 
The first one in the curly brackets is the difference between the perturbative analogue of the 
bilocal correlator and the continuum contribution into the "exact" one. This part cancels the 
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logarithmic non-analyticities in the perturbative diagram (eq. fl8|)) corresponding to the leading 
twist-2 operator in the OPE of pion currents. The tower of the next-to-leading non-analyticities 
can be subtracted in a similar way by accounting for the twist-4 operator. The last term displays 
the physical contribution to the correlation function that possesses the correct behaviour in the 
"momentum transferred Q 2 " . Requiring that in the limit of large Q 2 the expression in the braces 
should be zero, we come to the local duality relation for the overlaps 

(m M ) }j }j -( 1) 87r 2 2 ./-i r(2 j + 4)- W 
The sum of eqs. (||) and ( p5|) is an analytical function in Q 2 as all singularities are replaced 
by the combination Q 2 + Oj which is safe in the limit Q 2 —>■ 0. Due to the presence of the 
non-analyticities in each moment of the distribution function, we need an infinite number of 
parameters to be found from additional sum rules. Obviously, this is an impossible task. Safely, 
for the problem at hand, this part vanishes in the forward limit and the sum rule is dominated 
by the contact terms. 

The result of eq. (^), as concerns the Q 2 -independent part, can be seen from the WI. The 
contact term contribution contained in the bilocal correlator is effectively transformed into the 
power correction due to vector condensate which arises together with eq. (|8|) from the two-point 
correlation function in the WI (first two terms of eq. (|T2|)). The latter was investigated in 
connection with the pion wave function in the same framework J7|. However, it is not so for the 
most important bilocal part. In this respect the WI is useless as it transforms the bilocals which 
can be reduced to the condensates not accompanied by the strong coupling constant. 

The dominant contribution comes from the bilocal correlator convoluted with a three-propa- 
gator coefficient function (see fig. 1(e)) which looks like 

Wgl(plplQ 2 ,\) = /rf 4 xe^ x C| 2) (x)W|(g,x,An), (27) 

(the mirror conjugated contribution can trivially be added) where = Ax a + Bn a + Cpi a 
and we will not specify the coefficients in this decomposition because of their complexity. All 
nonperturbative information is accounted for in the correlator 

W^(g,x,An)=*| C /V 9y (0|T|w(0)$[0,x] M (x), o[y + ~n,y - ±n\\ |0). (28) 

In order to extract the contact term we make the following decomposition: 

W*{q, x, An) = + q^ 2) + An M Pf 3) (29) 

and convolute this expression with the vector q^. Performing the same steps as in the derivation 
of the Ward identity ( |1"2"D we find 



(qxjVfa = i(qx) J dre iT(qx) ( 



A A 

— , x 

2' 2 
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"4 / 1 1 ^/ rf V 9y (0|Tk(0)$[0,x]n(x),^(y+ \,y- ±n,A\ |0) + (A - -A) + Q 2 Pf 2) . (30) 



The last term in the second line vanishes in the limit of zero Q 2 , that manifests the absence of the 
massless particles in the corresponding channels. Within the accuracy we are limited to, we are 
left with the first term only because the second one contains an extra power of gG^, and thus 
the corresponding OPE starts from the higher orders in the coupling constant and the dimension 
of the operators. Performing the integration by using the method outlined at the beginning of 
the paper, we obtain the following contribution to the structure function: 

32 , fi 



(2) 9 2 f 

W BL {M , A) = — ixa s (uu) / dxcos(\x) 
9 Jo 



x / dyl dtl'dCfsl V 7 mAJ §M 2 W -f - , . (31) 



X J \C J \ x - v) N6C - xy£(\ 

Now, having accounted for additional terms in OPE, we can easily check that the normalization 
condition for the quark distribution in the pion is restored. 

7. For zero Q 2 the perturbative spectral density is concentrated on the line si = s%, so that 
there is no transition between the states with different masses. We collect all contributions and 
make the continuum subtraction that results in the substitution M 2 — > M 2 (l — exp(— so/M 2 )) 
in the perturbative term. We have found good stability of the distribution function with respect 
to the variation of the Borel parameter in the region 0.5 < M 2 < 0.8 for the standard value of 
the continuum threshold So = 0.7GeV 2 . The normalization point of the OPE is /i 2 ~ 0.5GeV 2 ; 
therefore, the function obtained can be regarded as an "input" quark distribution at this low 
energy scale. In fig. 2, we present the curves for the valence quark distribution in the pion for 
M 2 = 0.6GeV 2 : the solid and long-dashed lines correspond to the values of the average virtuality 
of vacuum quarks X 2 = O.QGeV 2 ( 7 " x = 0.15AGeV 2 ) and X 2 = OAGeV 2 (7- 1 = 0.087GW 2 ), 
respectively. In the large-x, region the corrections due to the quark condensate do not exceed 30% 
of the perturbative term. However, in the small-x region at x — 0.2 the ratio of the contact term 
to the main one comprises 50% for X 2 = 0.6Ge^ 2 and 70% for A 2 = OAGeV 2 . Below this point 
the nonperturbative contribution increases and reaches 100% at x = 0.13 for A 2 = OAGeV 2 (for 
A 2 = 0.6GeV 2 it still amounts 50%). So, for x as small as 0.2 we could not trust x-dependence 
of our result. Of course, there is no possibility to reproduce the correct x — > behaviour of 
the parton density in the present approach as it is determined by the exchanges of the Regge 
trajectories. 

Now we can comment on the contribution of the nonlocal gluon condensate to our sum rule. 
As can be easily seen from the WI some part of this contribution is concentrated in the two- 



point correlation function, which has been studied in ref. |21 |. Being numerically rather small, 
it contains terms not vanishing for distinguished from nonlocal quark condensates that 
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do not spoil the (1 — x)-behaviour as x — > 1, but only renormalize the slope. Therefore, the 
nonlocal gluon condensate limits the validity of the present approach from the large- X values. 
This conclusion is made discarding additional terms appearing from the three-point correlator in 
the WI which can somewhat change the situation. This problem, as well as a particular value of 
x in the large-x region, where the approach becomes invalid, deserves further investigation and 
only smallness of the gluon condensate contribution favours our decision to disregard it in the 
present study. 

Since our result is valid only in the limited region of Bjorken variable, we could not evolve 
it to the experimentally accessible energies. In fig. 2, we compare our calculation with the 



distribution obtained in the NJL model |22| at the same normalization point and find reasonable 
agreement between two approaches in a wide region of the momentum fraction. In fig. 3, the 
latter evolved up to Q 2 = 20GeV 2 (short-dashed curve) is compared with the presently available 
fits of experimental data |23| . It shows good agreement with the result of the analysis of Sutton, 
Martin, Roberts and Stirling (solid curve) |24| , which is consistent with all present Drell-Yan and 
prompt photon nN data. We also present the (long-dashed) curve due to Gliick, Reya and Vogt 
25| ; however, their result does not agree with E615 experiment [j23] which requires the valence 



distribution to be larger by 20%. Similarly enhanced distribution has been obtained in ref. [26 1. 
If the GRV curve is renormalized within a factor of 1.2 — 1.3 in the central region, there will be 
no disagreement between the different analyses. 

In conclusion, we have calculated the pionic parton density at low momentum scale in QCD 
sum rules with nonlocal condensates. It is shown that the parton sum rule is fulfilled only after 
the bilocal power corrections are accounted for. We have found good agreement with the u-quark 
distribution function computed in the NJL model which when evolved up to the experimental 
scales is well comparable with data. 
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Figure 1: Diagrams contributing to the operator product expansion of the correlation function 
(4): the first line display ordinary power corrections, while the second one — contribution due to 
the bilocal correlators. 




Figure. 2: Quark distribution in the pion at the low energy scale /i 2 ~ 0.5GeV 2 calculated 
from the QCD sum rule for different values of the average virtuality of vacuum quarks: solid and 
long-dashed curves correspond to A 2 = O.QGeV 2 and A 2 = OAGeV 2 , respectively. Short-dashed 
curve is the -u-quark density found in the NJL model [22]. 




Figure. 3: The experimental fits of the valence w-quark distribution in the pion at Q 2 = 
20GeV 2 : SMRS result [24] is depicted by solid curve, GRV analysis [25] is shown by long-dashed 
one. Short-dashed curve is the -u-quark density calculated in the NJL model [22] evolved up to 
Q 2 = 20GeV 2 . 



